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Abstract 

We present results from numerical evolutions of single black holes distorted 
by axisymmetric, but equatorially asymmetric, gravitational (Brill) waves. 
Net radiated energies, apparent horizon embeddings, and recoil velocities are 
shown for a range of Brill wave parameters, including both even and odd parity 
distortions of Schwarzschild black holes. We find that a wave packet initially 
concentrated on the black hole throat, a likely model also for highly asymmet- 
ric stellar collapse and late stage binary mergers, can generate a maximum 
recoil velocity of about 150 (23) km/sec for even (odd) parity perturbations, 
significantly less than that required to eject black holes from galactic cores. 



Typeset using REVTeX 



I. INTRODUCTION 



Many popular models of active galactic nuclei, quasars and even archetypical galaxies 
rely on the relativistic influence of black holes on the surrounding environment to provide 
power sources for observed spectral emissions and inferred motions of gaseous or stellar 
material. However, due to the strong gravitational effects of black holes, their role in the 
evolution of galactic cores and quasars is uncertain since numerical relativity computations 
are needed to perform detailed investigations of the near-field regime. In particular, grav- 
itational waves generated from sufficiently asymmetric systems (such as collapsing stellar 
cores and coalescing black holes) can carry a nonzero linear momentum component and 
impart a recoil velocity to the emitting objects, dominantly from the interplay between 
mass-quadrupole and mass-octupole or mass quadrupole and current-quadrupole contribu- 
tions [0,0. These velocities would be astrophysically significant if they were large enough 
to eject the emitting objects from the center of the host galaxy and send them hurtling 
through intergalactic space. Because the efficiency of momentum radiation emission is not 
known precisely, the dynamics and stability of systems containing black hole engines remain 
important but unresolved issues. If radiation reaction effects are significant, they may have 
considerable observable consequences for astrophysics and cosmology, including the redis- 
tribution and depletion of black holes from host galaxies, the disruption of active galactic 
core energetics, the introduction of black holes and stellar material into the intergalactic 
medium, and the general formation and structure attributes of galaxies. 

Although approximation studies of radiation recoil have been performed for more than 
two decades now, results from these calculations based on quasi-Newtonian and relativistic 
perturbation formalisms present an uncertain picture due to their incomplete treat- 
ment. More recently, Anninos and Brandt have numerically computed the recoil effect 
from fully general relativistic head-on collisions of two unequal mass black holes with time- 
symmetric initial data, and have shown that recoil velocities are of order 10 - 20 km/sec for 
black holes with moderately large initial separations ( ^ lOM, where M is the mass of the 
larger black hole). Their results are in rough agreement with, and generally confirm, various 
estimates from perturbation calculations. 

Here we continue to explore the radiation reaction process by computing the energies 
and recoil velocities from single black holes distorted by axisymmetric gravitational (Brill 
| |10| |) waves. We extend previous Brill wave + black hole investigations [|n],|T2| by relax- 
ing the equatorial mirror symmetry imposed in earlier work, thus allowing for mixtures of 
consecutive (even/odd) multipole contributions to the emitted radiation. In addition to 
investigating radiation reactions in this new class of physical systems, our implementation 
of Brill waves allows for very highly distorted black holes which can be thought of as mod- 
els for the late stage behavior of binary coalescing black holes. In fact, we have been able 
to simulate single black hole distortions ||T^JT3[| , as characterized by the ratio of polar to 
equatorial circumferences of apparent horizons, that are significantly greater than what we 
have observed in the merged state of two colliding black holes [|T3|-|15|] . Although the horizon 
distortion is not the only factor infiuencing recoil efficiency, one might nevertheless expect 
to obtain some idea, or perhaps even an absolute upper limit, of the recoil magnitude during 
the late stages of binary interactions by investigating strongly distorted single black hole 
systems. 



We generalize the prescription developed in references |]TT],|T2[ (and summarized for conve- 
nience in §|lD to specify equatorially asymmetric initial data and to parameterize Brill wave 
perturbations of Schwarzschild black holes by the amplitude, shape, location, and spectral 
mixture of the even and odd £-modes. Results from numerical evolutions are presented in 
§ [rn| , where we show embeddings of the black hole apparent horizons, energies emitted in the 
most dominant quasinormal modes of the final black hole, and recoil velocities arising from 
the mixing of consecutive radiative modes. The computations are carried out for both even 
and odd parity distortions of black holes, and over a wide range of wave strengths, initial 
placements, and mode distributions. We conclude in §ffV. 



II. INITIAL DATA 



For even parity distortions, we utilize the conformally fiat approach of Bowen & York 
to solve the initial value problem in axisymmetry and write the spatial 3-metric at the 
initial time as 
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where ^(77,6*) is the conformal factor, q{rj,6) is a function subject to certain constraints 
in its form but is otherwise freely specifiable, ^0(^5 ^) is chosen so that the Kerr metric is 
recovered if g = and the appropriate extrinsic curvature is specified, 77 is a logarithmic radial 
coordinate centered on the black hole throat, and {9, (p) are the usual angular coordinates. 
The more general 3-metric (applicable to both even and odd parity perturbations) is of the 
form 
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with F = (7^ 0) for the even (odd) parity cases. The remaining metric components (7^6) 
and ■jri(j>) are set to zero by the gauge freedom in choosing the shift vector. 

As in Ref. the somewhat arbitrary function q is restricted by symmetry conditions 
on the throat and axis, and fall-off rates at large radii [[l^]. The function q is constructed 
to have an inversion symmetric Gaussian part given by 



Qof{0) {e 



+ e 



where 
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and gi = or go for perturbations of the stationary Kerr solution or the Bowen & York 
spacetime |T6| respectively. With this form, the Brill waves are characterized by their 
amplitude Qq, width a, center coordinate location 770, and their angular dependence f{6). 
This allows a convenient way to parameterize the strength, shape and placement of the 
waves, and to easily tune the wave data for a broad range of spectral mode mixtures. 



In previous work, we had considered the case / = sin*^ 9 which possesses mirror symmetry 
across the equator. However, this data does not allow for the emitted gravitational waves to 
carry any linear momentum, as it excludes odd multipole components. Here, we relax the 
constraint of equatorial symmetry and consider 

/(n,e,^) = (l-e + ecos^)sin"^, (5) 

for which q in obeys the isometry conditions {j] —r], 9 —9, and 9 27i — 9). The 
form of (ID also has the necessary property /(O) = /(vr) = 0, and regulates the even and 
odd mode power distributions through the parameters n and ^. The parameter ^ determines 
the asymmetry of the wave, and the relative excitation of the odd and even numbered, even 
parity i modes. When ^ = and n = 2, £ = 2 is the dominant mode; when ^ = 1 and n = 2, 
£ = 3 is the dominant mode. For some intermediate value of ^ there will be a roughly even 
distribution of energy between the i = 2 and i = 3 modes, and at this value the gravity 
waves will produce a maximum recoil velocity on the black hole, as demonstrated in §|T|. 
The initial value problem for the even parity cases is completed by solving the Hamiltonian 
constraint for the conformal factor in the metric (|l|) with the specified free data. We 
also impose time symmetry, hence the extrinsic curvature is set to zero and the momentum 
constraint is trivially satisfied. This implies that the Brill wave packet is a combination of 
ingoing and outgoing radiation. 

For odd parity distortions, the free data is specified in the only nontrivial momentum 
constraint equation arising from the component in "time-rotation" symmetry and maximal 
slicing {tr K = 0). Defining the initial extrinsic curvature as 

He sin^ 9 ' 

^i^sin^ , (6) 

HEsm'^9 Hp sin 9 

the momentum constraint reduces to 

d^iHE sin^ 9) + deiHp sin^ 9) = 0, (7) 

and is satisfied by 

He = W) + Uv) [4COS0 U9) + sin^ dehm , (») 
HF = -dj2{v)s^r,^0 fsi9), (9) 

where /i, /2 and /s are arbitrary functions with the following symmetries: fi{9) = fi{—9) = 
/i(27r - 9), h{9) = fsi-9) = /3(27r - 9), = ^(-r/), and ^ ^ as r/ ^ oo so 

the spacetime asymptotically approaches the Kerr solution. To construct a model for odd 
parity waves analogous to the even parity case described above, we choose the following free 
functions: 

/i = 0, (10) 

/2 = Qo(e-^++e-^-), (11) 
/3 = (l-e + ecos^)sin"^. (12) 

Since He falls off sufficiently rapidly at large radii, any spacetime constructed using this 
conformal extrinsic curvature with /i = will have zero angular momentum. The Hamilto- 
nian constraint is then solved for ^ given the above extrinsic curvature and a conformally 
flat 3-metric with g = go = in (0). 
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III. RESULTS 



In this section we present results from several dozen numerical calculations of both even 
and odd parity axisymmetric distortions of single Schwarzschild black holes using the initial 
data parameterization described in §|T|. The results are presented as functions of various key 
Brill wave parameters, and their effect on radiation recoil is evaluated. In most cases we 
have used a numerical grid with 300x65 (radial x angular) zones to cover radial distances out 
to several hundred Madm-, where Madm is the ADM mass of the spacetime, and to include 
the entire polar domain < ^ < tt. However, we have also confirmed that the results are 
robust and relatively unchanged at different grid resolutions. Our simulations utilize the 
maximal slicing condition {K = K = and are generally run to 50 - TOM^^i^, which is 
more than enough time to extract the radiation content. The radiated wave energies and 
recoil velocities are computed from the energy-momentum flux across a spherical shell of 
radius I^Madm from the center of the black hole throat. 



A. Even Parity 

First we consider the effects of varying three independent parameters Qq, tjq and ^, 
corresponding to the amplitude, peak location, and dominant mode of the Brill waves, on 
the radiation reaction and dynamical evolution of even parity distortions of black holes. The 
remaining free initial data parameters described in §|I| have been held fixed: a = 1 for unit 
width wave profiles, qo = since we do not consider rotating black holes, and n = 2 to allow 
maximum grid resolution over the angular variations. 

From the equatorially symmetric examples in Refs. |TT]-|T^ , it is known that increasing the 
amplitude parameter Qq increases the strong field coupling of the Brill wave and black hole, 
substantially distorts the spacetime from spherical symmetry, and emits a greater fraction of 
the ADM mass in the form of gravitational radiation. To demonstrate the degree by which 
a black hole is distorted from sphericity, we first look at the geometric characteristics of the 
spacetime, namely the apparent horizon since it can easily be found in the spacelike slices. 
The horizon shape parameters and fiat space embeddings are evaluated for the case r/o = 
in which the Brill wave is placed directly on the black hole throat for maximum horizon 
distortion. We use a Newton-Raphson [|17| , p!8| procedure to solve the nonlinear equation 
defining the trapped surface conditions (zero expansion of outgoing null normals to the 
2-surface). The geometric properties are extracted from the two-dimensional sub-metric 
induced on the horizon surface 

d9^ + Ddcp"^ + 2Fd9d^^ , (13) 

where h{6) is the radial coordinate defining the horizon. Visual representations of the horizon 
are achieved by embedding the 2-surface given by (p!3| ) in a higher three-dimensional fiat 
space. Introducing a new coordinate z on a fiat 3-metric, the 2-metric of the horizon surface 
is identified as 
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where B' = + A{dh/d6)'^ + F'^/D) and D' = '^'^D are the metric components of the 

horizon surface transformed to a diagonal form. Solving for the coordinate z gives 



z = J de^B' - {de^')\ (15) 

which is integrated numerically to obtain the embedding functions z{0) and p{6) = y/D', 
although an embedding is not in general guaranteed to exist. 

In figure we show embeddings of the horizon 2-surface in the initial data (where 
distortions are greatest with 770 = 0) for a highly perturbed Qo = 0.9 case with different 
mode parameters ^. Embeddings of the more prolate odd mode distortions, i.e. ^ > 0.7, 
are undefined (and therefore not displayed) for the more negative values of z due to the 
radical in equation (|T5|) which becomes negative. Figure ^ shows the equivalent embeddings 
for ^ = 0.5 as a function of wave amplitude Qq. These embedded distortions eventually 
damp out in an oscillatory fashion over time as the horizon evolves towards sphericity after 
the dynamic component has either been absorbed by the black hole, or has propagated to 
asymptotic infinity in the form of gravitational waves. Together, the embedding diagrams 
indicate that horizon distortions are roughly spherical for small amplitude perturbations 
and become generally more prolate as Qo is increased. The shape of distortions also varies 
with the mode parameter ^, which regulates changes from equatorially symmetric even mode 
behavior, to predominantly asymmetric odd mode configurations for the larger values of ^. 
The ratio of polar to equatorial circumferences of the horizon surfaces, used in our previous 
work as a measure of distortion, is not especially informative regarding the magnitude of 
radiation recoil. Indeed, the purely even symmetric cases generally give rise to greater 
distortions, but no recoil which is a function of the relative mixture of even and odd modes, 
as well as the perturbation amplitude. Assuming a simple definition of radial distortion 
in the embeddings as Rr = max(-\/p^~+^)/min(v^p^~-i-^), the displayed distortions range 
from Rr = 1.2 for (Qo = 0.1, ^ = 0.9), to 7.5 for (Qo = 1-2, ^ = 0.1). 

The transition from even to odd mode behavior observed in the near-field horizon em- 
bedding diagrams of figure |I|, is also mirrored by the mode distribution in the far-field 
radiation zone. Figures ^ and ^ show the energy (normalized to the ADM mass of the 
spacetime) radiated in the most dominant i = 2 and 3 modes for r^o = as a function of ^ 
and Qq. As the mode parameter is increased in figure ^ the energy distribution dominance 
changes from even to odd, consistent with the horizon embeddings in figure ^ Figure ^ 
indicates that the total radiated energies asymptotically approach constant maximal values 
for each of the mode parameters, and that the range of parameters we have investigated are 
reasonably representative of the most efficient radiators of gravitational energy. We have 
restricted current studies to wave amplitudes Qo < 1.2, since the numerical results are less 
reliable for larger amplitudes, especially at late times and in the ability to resolve both the 
i = 2 and 3 modes in the extreme odd or even i-mode dominated evolutions. 

The mixing of adjacent multipole modes gives rise to a non- vanishing flux of linear 
momentum along the z-axis which can be evaluated from products of consecutive Zerilli 
wave functions IHI 
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where the Zerilh functions ijjg are normahzed such that the total radiated energy in each 
mode is given by 



^'=32^ J "^^^^'^ • ^^^^ 

For numerically practical purposes, we compute only the most significant (2,3) and (3,4) 
contributions. In general we find that the higher order terms in the series (|16|) play an in- 
creasingly greater role as ^ is increased and as the distortions become dominantly odd func- 
tions. For an interesting large amplitude case (Qo = 0.9) the momentum ratio -P(2,3) /-P(3,4) 
varies from roughly 50 to 0.2 for ^ = 0.1 and 0.9 respectively, with even greater ratios for 
the smaller amplitude cases. However, as we show below, the greatest recoil velocities arise 
for roughly equal mixtures of even and odd perturbations (i.e., C, = 0.5), and for these cases 
the (2,3) contribution exceeds the (3,4) by at least a factor of 50 in all cases we have studied. 
The results presented for the radiated momentum are derived by summing both the (2,3) 
and (3,4) contributions. To confirm the degree to which these two dominant pairs are com- 
plete, and to independently check our calculations, we have also evaluated the momentum 
fiux from the Landau-Lifshitz pseudotensor ||T9|-pl| 
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cosOdn, (18) 



which is valid for both even {F = 0) and odd {B = D = 0) parity perturbations. We find 
excellent agreement, typically better than 10%, between the two calculations. 

As a result of the momentum emission from the interplay between even and odd £-modes, 
the final black hole will acquire a recoil velocity 

(19) 
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that is opposite in direction to the momentum fiux of the waves. We show these velocities in 
figures 1^ and ^ in physical units of kilometers per second, and as a function of the parameters 
^, rjo and Qq. Figure ^ confirms that maximum recoil occurs at ^ = 0.5. The asymptotic 
fiatness of the curves at large values of Qo in figure ^ suggests maximum velocities of about 
150, 200, and 500 km/sec for the most strongly distorted cases that we are able to compute 
numerically for r]Q = 0, 0.25 and 0.5 respectively. However, we are unable to reliably evolve 
greater amplitudes for rjo = 0.25 or 0.5 and establish a precise turnover velocity since our 
code breaks at these large amplitudes, though the curves already show signs of fiattening by 
Qq = 1.2. Hence the quoted values in these two cases are approximate extrapolations. Also 
notice that greater recoil velocities result when the Brill waves are placed at larger radii, 
beyond the perturbation potential barrier (at r > 3M, or equivalently rj > 1.8). In these 
cases, the ingoing waves excite the ringing modes more strongly as they cross the barrier 
and emit a greater fiux of energy-momentum. On the other hand, it is also likely that a 
large fraction of the radiated fiux in the large tiq cases can be attributed to the outgoing 
wave component and general distortions of the global spacetime [0, and not to pure black 
hole ringing from localized collapse or impact scenarios. 



B. Odd Parity 



Following the general presentation of §|III A|, we present in this section recoil velocities 



from odd parity distortions of black holes as a function of mode parameter (^), initial peak 
location (r/o), and wave amplitude (Qq) in the initial data of §|I[ 

The momentum arising from consecutive £-mode interactions of odd parity waves takes 
a form analogous to ([16|), except the wave functions ip'^'^'^ (replacing the even parity ipe) are 
extracted from the 'je^ metric component by 



re 



2)\ 



\ (^ + 2)! 



J d^r^F^dl - cot dde)Y(,QdVl. (20) 



In this form, the wave functions are related to the Regge- Wheeler perturbation variable 
h2 



with normalization 
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An important difference in the evolutions of odd (versus even) parity data is that it 
is necessary to keep track of the momentum contributions from a greater number of i- 
mode pairs since there is no clear dominance by the lowest order terms. The individual 
contributions can jump from positive to negative values of momentum with roughly the 
same amplitude over consecutive pairs. For example, the better behaved r/o = small 
amplitude cases require at least three pairs in the series to converge at the 10% level when 
compared with subsequent additions and to the Landau-Lifshitz formula (|T^). Also, the 
odd parity initial value problem and the dynamical evolutions over time can both generate 
significant even parity signals, contributing up to a few percent of the net recoil velocity 
in the cases we have investigated. Hence all results in this section are derived from the 
Landau-Lifshitz pseudotensor, thus accounting for all the modes in evaluating the velocities 
(though we have compared and confirmed the consistency of both methods for small and 
large amplitude runs). We plot in Figures |^ and ^ the velocities as a function of ^ and 
Qo respectively. Figure ^ indicates that maximum recoil is achieved for ^ 0.7. Figure |^ 
shows the maximum velocity for ^ = 0.7 as Qo is varied over a numerically robust range 
of amplitudes for three different initial wave positions rjo = 0, 0.5 and 1. The evolutions 
generate maximum velocities of 23, 52 and 430 km/sec, with increasingly greater velocities 
for Brill waves initially concentrated at greater distances from the black holes. 

In comparing these results with the even parity cases (say tjq = 0), one should not 
necessarily conclude that odd parity radiation is less effective in producing radiation recoil 
for any intrinsic reason. Much has to do with the manner in which the data was constructed. 
For example, the even parity data distorts the spatial metric, while the odd parity data uses 
a conformally flat metric. One could, in principle, produce even parity distortions of the 
spacetime through the extrinsic curvature while maintaining a fiat metric. Likewise, one 



could add odd parity data to the metric itself rather than the extrinsic curvature. This would 
make the procedures more similar (and possibly the radiation energies as well). Furthermore, 
there are many ways to construct initial data for both types of radiation and it is not feasible 
to study them all. Rather, our results represent the maximal effects of a certain class of 
black hole distortions. 

IV. CONCLUSION 

We have carried out a systematic study of single black holes distorted by strong-field 
axisymmetric Brill waves in an effort to quantify the astrophysical significance of the "rocket" 
effect imparted to the final black hole from the momentum carried by gravitational radiation 
in the system. This work compliments our previous studies of the head-on collision of 
two unequal mass black holes where we found recoil velocities up to 10-20 km/sec. 
However, it is likely that coalescing binary black holes with arbitrary physical parameters 
(i.e., impact parameters, masses, and spins) may generate greater recoil velocities, so we 
have focused these current studies to deduce the maximum recoil expected from highly 
asymmetrical configurations. The Brill wave + black hole systems we have studied allow 
a parameterization of the wave strengths, widths, locations, and shapes of the perturbing 
sources such that we are able to systematically explore the role of various parameters in 
fully nonlinear numerical calculations of strongly distorted black hole spacetimes. With this 
approach, we are able to generate greater distortions and wider spectral energy distributions 
of black holes than observed in our simulations of colliding binary systems. We thus also 
consider our current results as reasonable maximum estimates of radiation recoil in single 
or late-stage binary black hole systems (although a more precise comparison between single 
and binary evolutions must also account for any residual radiation content in the initial data 
of the respective systems). 

For the most highly distorted spacetimes, we find maximum recoil velocities in excess of 
400 km/sec for both even and odd parity data with Brill waves initially centered at large 
distances from the black hole throat, e.g., rjo = 0.5 (1.0) for even (odd) parity perturbations. 
Our results exhibit a strong dependence on the initial placement of the Brill waves, as well 
as their amplitude and spectral composition. Of all these effects, we are less certain of 
the role which the initial wave placement t]o plays in generating a true maximum value, 
since for the numerically difficult combination of large separations and amplitudes, our code 
eventually breaks down. However, we expect for radiation clumps located further from the 
black hole, that a substantial fraction of the (outgoing component of the) Brill waves escapes 
to infinity since the perturbations are applied essentially to the spacetime surrounding the 
black hole, and not directly on the throat. Hence we expect that the bulk of emitted energy- 
momentum flux can be attributed to the initial wave configuration for large rjo, as opposed 
to any intrinsic ringing of the black hole associated with localized source dynamics such as 
ingoing wave collisions, collapsing stellar cores, or coalescing binaries. On the other hand, 
an ingoing wave located outside the potential barrier can scatter off and impart a much 
greater momentum to the hole. We were, however, unable to distinguish secondary wave 
pulses in our numerical data corresponding to reflected waves. 

In addition to reducing the initial outgoing Brill wave content, it is also likely that the 
rjQ = cases represent more appropriate late stage recoil models for black hole binary sys- 



terns. In these cases, our results of 150 and 23 km/sec for even and odd parity perturbations 
are in general agreement with the bound Vr < 300 km/sec derived by Bekenstein in his 
quasi-Newtonian considerations of the interaction between quadrupole and octupole terms 
in non-spherical stellar core collapse to black holes. Furthermore, the odd parity recoil in our 
calculations is remarkably similar to the 25 km/sec found by Moncrief [^] for non-spherical 
models of black hole formation. Our even parity results are approximately a factor of two 
times larger than the quasi-Newtonian calculations {vr ~ 67 km/sec) of binary systems in 
Keplerian orbits by Fitchett p. However, considering the ambiguity in choosing the final 
prior-to-plunging orbit and in extrapolating the perturbation calculations to the equal mass 
limit, our results are in fairly good agreement with the predictions of Fitchett and Detweiler 
0] who extended Fitchett's earlier work to perturbation theory and computed a maximum 
velocity of about 120 km/sec for the merging of two black holes from the last stable circular 
orbit. We are also in agreement with Nakamura, Oohara and Kojima who estimate a 
maximum velocity of about 240 km/sec from numerical perturbative calculations of test 
bodies plunging into black holes from infinity with arbitrary orbital angular momentum. 

In comparison, the escape velocity from galactic structures can vary from about several 
hundred km/sec for spiral galaxies such as the Milky Way, to about a thousand km/sec for 
the more massive giant ellipticals such as M87. Our results, however, suggest that black 
holes which may be located in the centers of galaxies and which undergo highly asymmetric 
evolutions (including strong field distortions and binary mergers) are relatively stable entities 
and will not likely escape from the host galaxy, assuming that the "on the throat" numerical 
calculations are reasonably representative models. Although we have established that the 
recoil effect is not generally large enough to be considered astrophysically significant, this 
does not, however, rule out the possibility of black hole ejections from galactic disks far from 
the core and in the direction of galactic rotation, nor the possibility of black hole ejections 
from globular cluster systems in galactic halos. Black holes can more readily escape from 
these systems to wander through the galaxy or even intergalactic space. 
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FIGURES 




FIG. 1. Embeddings in 3-dimensional flat space of tlie apparent liorizon 2-surfaces at tlie 
initial time for even parity perturbations. Each curve represents the p-z embedding for the same 
strongly distorted case with Brill wave amplitude Qq = 0.9 and wave location r]Q = 0, but with 
different mode parameters ^ that regulate the relative mixture of even and odd mode spectral 
distributions. Distortions range from i?^ = 1.9 for the ^ = 0.9 case to Rr = 4.9 for ^ = 0.1, where 
Rr = max(-\/ p2 -|- 2;2)/niin(\/ + z^) is the ratio of the maximum to minimum distances of the 
embedding curves to the origin at /? = z = 0. 
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FIG. 2. As figure |T|, except here each curve represents the embedding for the same mode pa- 
rameter ^ = 0.5 which generates the greatest recoil effect for a specified amplitude Qq. Distortions 
range from Rr = 1.2 for Qo = 0.1 to Rr = 2.7 for Qo = 0.9. 
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FIG. 3. Radiated energies (normalized to the ADM mass, Madm) in the most dominant 
£ = 2 and 3 even parity modes as a function of the Brill wave mode parameter ^ for 7/0 = and 
select values of Qq. The wave extractions are performed on spherical shells of radius I^Madm and 
centered on the black hole throat. 




FIG. 4. As figure except the energies are plotted as a function of Brill wave amplitude for 
select values of The upper bound on the wave amplitudes considered here is dictated by the 
accuracy with which the weaker of the £ = 2 or 3 modes can be resolved. 
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FIG. 5. Recoil velocities for even parity distortions as a function of mode parameter ^ for 
ijQ = and different wave amplitudes Qq. The velocities are computed from the radiative momen- 
tum flux carried by gravitational waves crossing spherical shells of radius 15Madm from the black 
hole throat. 




FIG. 6. Recoil velocities for even parity perturbations as a function of Brill wave amplitude 
Qo for ^ = 0.5 and different initial wave locations rjQ. The trend toward asymptotic flatness in the 
displayed curves suggests maximum recoil velocities of about 150, 200 and 500 km/sec, for r/o = 0, 
0.25 and 0.5 respectively. 
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FIG. 7. Recoil velocities for odd parity distortions as a function of mode parameter ^ for 
rjo = and different wave amplitudes Qq. 
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FIG. 8. Recoil velocities for odd parity perturbations as a function of Brill wave amplitude 
Qo for ^ = 0.7 and different initial wave locations 770. Maximum recoil velocities of about 23, 52 
and 430 km/sec are found for r/o = 0, 0.5 and 1 respectively. 



